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ABSTRACT 

Aims. We examine the physical conditions of ionised gas in the central kiloparsec(s) of nearby Seyfert galaxies. Our aim is a detailed 
two-dimensional mapping of optical spectroscopic properties, searching for features in the emission structure and in gas kinematics 
which would be common to all the active galactic nuclei (AGN) or their class, and which would help us to understand the origin and 
nature of the circumnuclear gas, and its role in the AGN. 

Methods. We map narrow-line regions (NLRs) of 1 1 nearby Seyfert 2 galaxies with the optical integral-field spectrograph OASIS 
mounted at CFHT. We model emission-line profiles of 5 forbidden-line doublets and 2 Balmer lines, correcting for the underlying 
stellar absorption by reconstructing stellar spectra with synthetic evolutionary stellar population models. 

Results. For each of the 1 1 targets, we present 2D maps of surface brightness in the observed emission lines, diagnostic line intensity 
ratios, gas kinematics (mean line-of-sight velocity and velocity dispersion), electron density, and interstellar reddening, and we plot 
spatially resolved spectral-diagnostic diagrams. The stellar data are represented by maps of mean line-of-sight (LOS) velocities 
and of the relative mass fractions of the young stellar populations. The gas velocity fields in 80% of the sample exhibit twisted S- 
shaped isovelocity contours, which are signatures of non-circular orbits and indicate non-axisymmetric gravitational potentials, gas 
motions out of the galactic plane, or possible outflows and inflows. Based on the kinematic measurements, we identified a possible 
nuclear ring or radial gas flow in NGC 262 (Mrk 348), not reported before. Eight of the eleven observed objects exhibit strongly 
asymmetric or multi-component emission-line profiles, in most cases confined to an elongated region passing through the galactic 
centre, perpendicular to the major axis of emission. 

Key words, galaxies: active; galaxies: Seyfert; galaxies: kinematics and dynamics; methods: observational; methods: data analysis; 
techniques: spectroscopic 



1. Introduction 

In the search for the fuelling mechanism of active galactic nuclei 
(AGN) and its trigger, no significant differences have been found 
between the host galaxies of AGN and their quiescent counter- 
parts. Among others, no excess of galactic bars has been con- 
firmed in active galaxies (Martini et al.|2003 |Laurikainen et al. 
|2004) . On the other hand, recent studies of galactic centres have 
reported an excess of infrared isophotal twists ( Hunt & Malkan 
and of twisted isovelocity contours in gas (Dumas et al. 
which may imply presence of non-axisymmetric poten- 
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tials ( Schinnerer et al.|2 000). Focus on gas kinematics in the nu- 
clear regions thus seems to be the key to understanding the inter- 
play between the host galaxy and the active nucleus, even though 
interpretation of the phenomena on very different spatial scales 
needs to be done with extreme care. Detailed spatially resolved 
spectroscopic observations of the nuclear regions in large sam- 
ples of galaxies, followed by proper modelling are hence neces- 
sary. 

The Narrow-Line Regions (NLR) represent low-density gas 
ionised by the central AGN (e.g., [Morse et al.|1996[ Evans et al. 
|1999| l and extending beyond the central dusty torus. The clas- 



sical NLRs lie within 10 - 100 pc scales and their kinemat- 
ics is often perturbed by energy injection from the AGN while 
gas motions in the so-called extended NLRs (ENLRs), reach- 
ing to kpc and 10 kpc scales, are thought to be predominantly 
gravitational. The ENLRs are convenient observational targets 
for spatially resolved imaging and spectroscopy, especially in 
nearby Seyfert galaxies. They play a crucial role in attempts to 
understand the AGN structure and evolution, and in testing of 
the Unified Models of AGN ( [Antonucci 1993 ; Urry & Padovain 
J995]l. 

By completing narrow-band imaging and aperture spec- 
troscopy observations of large samples of Seyfert galaxies, ro- 
bust statistical analyses of NLR properties have been achieved 
(Mulchaey et al.||1996a[ |Ho et al.|[Tg97l [Schmitt et al.]|2003a| l. 
Similar surveys with modem spatially resolved spectroscopic 
techniques are necessary to obtain more complex information 
and achieve a qualitatively new level of conclusions. As an ex- 
ample, we mention the controversy in optical imaging studies by 
Mulchaey et al.| ( |1996al l; and |Schmitt et al.| ( [2003a) l, who tested 
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predictions of the Unified Model with respect to the NLR size 
and morphology, and obtained results that cast doubts on the 
simple model of ionisation of ambient gas by a cone of radia- 
tion. However, both studies are affected by the limitations typi- 
cal of pure imaging: these include the survey sensitivity or con- 
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tamination by other emission sources, the most important be- 
ing gas ionised by young, massive stars that have been found in 
centres of many Seyfer t galaxies (e.g., Schmitt et al. 1999; Cid] 
[Fernandes eFa l. 2004; Boisson et al.|2004| l. " ' 

Spatially resolved spectroscopy overcomes many of these 
difficulties: ratios of forbidden and permitted emission-line in- 
tensities at different locations help us to discriminate between 
the ionisation sources and to establish a boundary of the AGN 
ionised NLR, using ionisation models (see e.g. 



Bennert et al. 



2006a| for a discussion and application). Also, spatially resolved 



kinematics derived from the emission lines provides an impor- 
tant insight into the three-dimensional structure and dynamics 
of the central emitting regions. The wealth of information con- 
tained in such spectroscopic data thus allows us to address a 
number of questions related to the NLRs, which are important 
for reconstructing the global picture of AGN. An efficient way 
of acquiring simultaneously the spatially resolved optical or in- 
frared (IR) spectroscopic data is provided by the integral-field 
spectroscopy (IFS). Integral-field units (IFUs) are becoming a 
standard type of equipment at large telescopes, available with 
a wide variety of technical parameters, such as spectral range, 
spectral and spatial resolution, and field of view. 

We present our mapping of the central ~ kiloparsec(s) of 
1 1 nearby Seyfert 2 galaxies observed with the optical integral- 
field spectrograph OASIS mounted at the 3.6 m Canada-France- 
Hawaii Telescope (CFHT). Now mounted at William Herschel 
Telescope (WHT), OASIS is an IFU well-suited to high- 
resolution observations over a restricted field of view (FOV), for 
a wide range of spectral modes. Its use for detailed spectroscopic 
mapping of the NLRs and the central parts of the ENLRs is thus 
advantageous for understanding the complex spatial variations 
in the physical properties of the circumnuclear ionised gas, and 
the implications for the AGN model. 

For the observed emission regions, we plot maps of derived 
quantities such as the surface brightness in two Balmer lines 
and five doublets of forbidden emission lines, the computed line 
spectral-diagnostic ratios, the kinematic properties, the electron 
densities, and the interstellar reddening. We plot spatially re- 
solved spectral-diagnostic diagrams for each object where pos- 
sible. We provide details of the individual Seyfert 2 nuclei in 
Sect.[4j confronting our observations with results found in the 
literature. The complete sample is discussed in Sect.|5] where we 
draw conclusions about various structural and kinematic results. 



2. Observations and data reduction 

2.1. Observations 

The selected 11 targets (Table[T]i correspond to well-known 
Seyfert galaxies with previously published photometric and/or 
spectroscopic data. Their redshifts are in the range 0.006 < 
z < 0.051 , corresponding to distances of 17Mpc - 220 Mpc. 
Observations were performed in the years 2000 - 2002, with the 
lenslet-array OASIS spectrograph constructed at the Centre de 
Recherche Astrophysique de Lyon (CRAL), France, mounted 
at CFHT (Mauna Kea) in the F/8 Cassegrain focus. Two spec- 
tral domains were explored, with either high (0.27") or low 
(0.41") spatial sampling, and corresponding fields of view of 
10" x8" and 15" x 12", respectively. Each of the modes provided 
approximately a thousand spectra per object per spectral do- 
main. The two spectral modes covered the wavelength ranges of 
4760 A-5558 A with 2.15 A dispersion ("MRl" spectral mode). 



and 6210A - 7008A with 2.17A dispersion ("MR2" spectral 
mode). 

Table[2] specifies the spectral and spatial configurations used 
for each member of the entire observed set. An overview of all 
the configurations of OASIS is available at the dedicated web 
pages of the instrument and at the CFHT pages. Details of the 
instrumental concept are to be found in Bacon et al. ( 1995[ l, 
which describes the TIGER spectrograph, a direct predecessor 
of OASIS. 



2.2. Data reduction 



The data were reduced with ffie dedicated XOasis software 



( Pecontal-Rousset et al. 2004a b). The standard reduction pro- 



cedure includes bias and dark subtraction (CCD), extraction of 
the spectra using a fitted mask model, wavelength calibration, 
low-frequency flat-fielding, spectro-spatial flat-fielding, cosmic- 
ray removal, homogenisation of the spectral resolution over the 
field, sky subtraction, and flux calibration using observations of 
the photometric standard stars. Multiple exposures were merged 
and mosaiced, truncating the wavelength domain to a common 
range in the different fields, and combining the spectra (both sig- 
nal and noise) with optimal weights and renormalisation. The 
data cubes were resampled to a common spatial scale on a rect- 
angular grid (from the original array of hexagonal lenses) with 
either 0.25" or 0.4" spacing depending on the adopted spatial 
sampling. Orientation was defined to be north up, east to the left. 



3. Data analysis 

The results presented in this paper were obtained by a two-stage 
analysis procedure: stellar-population modelling and emission- 
line modelling. The OASIS survey that we present was designed 
preferentially for the study of the gas dynamics, and the signal- 
to-noise ratio (S/N) of the acquired data was therefore calculated 
for optimal results in the gas emission lines. A simplified stellar- 
population modelling, with the primary goal of correcting the 
Balmer lines of hydrogen for underlying stellar absorption, per- 
mitted us to investigate both the stellar kinematics and the age 
of stellar populations. 

Methods for the simultaneous fitting of stellar and emission 
gas spectra were recently developed by the SAURON consor- 
,|Sarzi et al.|2006[|McDermid et al.|200 6) and oth- 



tium (se e e.g., 

ers (e.g., [Sanchez et al.|2007 1. We briefly tested this software on 



our data and found that the low S/N in our stellar lines prevented 
us from judging reliably which method was more suitable for 
our data. In addition, the two-stage spectral analysis, with the 
independent treatment of the stellar and gaseous components, 
is advantageous for detailed studies of complex emission-line 
profiles: the stellar component removed, different emission-line 
models with an increasing degree of precision can gradually be 
applied, as needed. 

For each of the data-analysis steps, the large volumes of data 
require to a large extent automatic treatment, with a subsequent 
check of the results, and detailed modelling applied only to se- 
lected spectra. The current paper presents our first-order results, 
while detailed analysis and interpretation of selected objects is 
postponed to future papers dedicated to individual objects in the 
sample. 
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UDject 


Alternative 


KA 




rlUDDle 


Keasniit z 


Velocity cz 


r\ 
Ul 


Ang. scale 




name 


[h m s] 


[d m s] 


type 




[km s ^ 


1 

J 


[MpcJ 


Lpc/ J 


Mrk 34 


MCG +10-1 5-104 


10 '^4 08 6 


+60 01 S2 


Spiral 


O'SO'S 


15 140 


+ 90 


218 


956 


Mrk 622 

±V±±JV \J^^ 


T IGC 04229 


08 07 4 1 

\J\J \J 1 i- .\J 


+39 00 1 5 

\J\J 


SO 


0.0232 


6964 ± 


11 


99.6 


461 


Mrk 1066 


UGC 02456 


02 59 58.6 


+36 49 14 


(R)SB(s)0+ 


0.0120 


3605 ± 


22 


47.3 


224 


NGC 262 


Mrk 348 


00 48 47.1 


+31 57 25 


SA(s)0/a 


0.0150 


4507 ± 


4 


58.2 


274 


NGC 449 


Mrk 1 


01 16 07.2 


+33 05 22 


(R')S? 


0.0159 


4780 ± 


2 


62.3 


293 


NGC 2273 


Mrk 620 


06 50 08.6 


+60 50 45 


SB(r)a: 


0.0061 


1840 ± 


4 


25.8 


124 


NGC 2992 


MCG -02-25-014 


09 45 42.0 


-14 19 35 


Sa pec 


0.0077 


2311 ± 


4 


36.6 


175 


NGC 3081 


IC 2529 


09 59 29.5 


-22 49 35 


(Ri)SAB(r)0/a 


0.0080 


2391 ± 


3 


37.7 


180 


NGC 4388 


UGC 07520 


12 25 46.7 


+ 12 39 44 


SA(s)b: sp 


0.0084 


2524 ± 


1 


I6.7W 


81 w 


NGC 5728 


MCG -03-37-005 


14 42 23.9 


-17 15 11 


(Ri)SAB(r)a 


0.0094 


2804 ± 


20 


41.9 


199 


NGC 5929 


UGC 09851 


15 26 06.1 


+41 40 14 


Sab: pec 


0.0083 


2492 ± 


8 


35.7 


170 



Adopted from the NASA/IPAC Extragalactic Database (NED). The cosmological scales assume Ho = 73kms"' Mpc"', £1„ = 0.27, Qa = 0.73 . 
Dt. luminosity distance; Ang. scale: angular scale. (*) Values from |Yasuda et al.| ( |1997) , with NGC 4388 considered inside the Virgo cluster. 



Table 2. List of observed Seyfert galaxies - observational parameters. 



Object 


Observation 


Integration time 


Integration time 


Spatial sampling 


Spatial sampling 




date 


in MRl [s] 


in MR2 [s] 


MRl 


MR2 


Mrk 34 


14-15/03/2001 


1 X 2700 


2 X 2700 


HR 


LR 


Mrk 622 


24/1 1/2000 


2x 1800 


2 X 1800 


HR 


HR 


Mrk 1066 


22/1 1/2000 


3 X 1800 + 1 X 1200 


2 X 1800 


HR 


HR 


NGC 262 


25/1 1/2000 


2x 1800 


2 X 1800 


HR 


HR 


NGC 449 


24/1 1/2000 


3 fields x2x 1800 


3 fields x2 X 1800 


HR 


HR 


NGC 2273 


24/1 1/2000 


2x 1800 




LR 




NGC 2992 


16-17/03/2001 


1 X 1800 + 1 X 2700 


2 X 1800 


LR 


HR 


NGC 3081 


17/01/2002 


2 X 3600 




LR 




NGC 4388 


17/01/2002 


3 X 3300 




LR 




NGC 5728 


15/03/2001 


2x 1800 


1 X 1800 + 1 X 2700 


LR 


LR 


NGC 5929 


16/03/2001 


2 X 2700 




LR 





The abbreviations used: MRl : spectral mode of wavelength range (4760, 5558) A, dispersion 2.15 A/pixel, resolving power 1210; MR2: spectral 
mode of wavelength range (6210, 7008) A, dispersion 2.17 A/pixel, resolving power 1525; HR: high spatial resolution mode with sky sampling 
0.27", field of view 10.4" x 8.3"; LR: low spatial resolution mode, sky sampling 0.41", field of view 15" x 12". 



3.1. Stellar populations 

We modelled the stellar component of spectra using synthetic 
evolutionary models from the GALAXEV library of Bruzual &] 



Chariot ( 2003| l. Understanding the stellar populations in Seyfert 
galaxies has always been a challenging problem (see e.g., |Cid| 
Fernandes et al. 2005 ; Raimann et al.|2003 ), due to the presence 
of strong gas emission lines, which contaminate the stellar spec- 
tra. A fairly standard approach has been to correct emission lines 
for stellar absorption by assuming that the nuclear stellar popu- 
lations are identical to those in the outer parts of the host galaxy 
(e.g., |Ho et al.|1997[|Bennert et al . 2006a). However, evidence of 
young stars in the centres of Seyferts has been presented by e.g., 
Heckman et al. ( 1997 1; Gonzalez Delgado et al. ( 1998|l; Joguet 



et al. 



( |2001| l; [Raimann et al.| ( |2003| l; and |Cid Fernandes et al. 



( 2004[ i, and therefore the use of purely old stellar populations 
has been shown to be insufficient even for detailed gas-emission 
studies, because of the need to correct reliably for the effects of 
stellar absorption. 

Stellar-population modelling searches for a combination of 
template spectra from a theoretical or experimental library which 
reproduces best the observed galactic spectrum. As a typical in- 
verse problem, the stellar spectrum reconstruction has a number 



of free parameters (such as the ages and metallicities of the indi- 
vidual populations, and their mass fractions), which may not be 
possible to constrain to form a unique solution. Intrinsic degen- 
eracies exist between age and metallicity, and they are further 
enhanced by measurement uncertainties and noise, which wash 
away the differences between spectrally similar components, as 
demonstrated by Cid Fernandes et al.| ([2004). These authors 
showed for artificial data that even an accurately fitting template 
does not ensure a correct physical interpretation. Nevertheless, 
a most robust determination of population fractions is possible 
if larger stellar-age bins are used (defined coarsely as young 
< 100 Myr, intermediate 100 Myr - 1.4Gyr, and old > 1.4Gyr 
in their study). 

We adapted the stellar population study to our data quality. 
We were limited by the S/N in stellar absorption lines, the rela- 
tively short wavelength ranges with only a few prominent ab- 
sorption lines, and contamination by gas emission lines. The 
fitting of the stellar was only possible in one of the two ob- 
served spectral regions, i.e., the MRl domain (see Sect. 2.1 1, 
which contains the Mgib complex at about 5 170 A, and iron 
lines Fe 5015 A, Fe 5270 A. Other constraints to the models were 
provided by the absorption-line wings of H/3. No such features 
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were detected in the MR2 domain, and therefore the stellar com- 
ponent of spectra in this domain was inferred by extrapolating 
the models fitted to MRl. 

From the GALAXEV models, we selected those that rep- 
resent evolution of Simple Stellar Populations (SSPs) with 
Salpeter initial mass function (IMF), by applying Padova evo- 
lutionary tracks (e.g., Girardi et al. 2000). We based our work 
on the findings of |Tremonti et al. (2004 ), and |C id Femandes| 
|et al.| ( |2004) l, who tested different sets of SSP models for recon- 
structing spectra of active and non-active galaxies. They defined 
a spectral base of 30 SSPs covering ten ages and three metal- 
licities, which played a role of independent functions, and any 
other SSP was found to be their linear combination. The tests 
that we performed with this "Tremonti set" on our data revealed 
that the number of free parameters introduced by the theoreti- 
cal combinations within the 30 SSPs was too large for our data 
quality, and could not be sufficiently constrained to produce a 
unique model, resulting in significant degeneracy with a difficult 
interpretation. We therefore searched for a simpler representa- 
tion, and used models combining two SSPs of markedly differ- 
ent ages, as applied before e.g., by Canalize et al. (2007)). The 
even more basic option of single-population models proved to be 
inadequate - the models failed to reproduce the observed spec- 
tral features correctly. The two SSP models were selected to be 
of solar metallicity, and of ages of 11 Gyr and lOOMyr, which 
were found to contribute the most even in the models to our data 
produced by the entire Tremonti set, and were able to describe 
many different types of spectra across the whole sample. The fi- 
nal fits differ from those obtained with the Tremonti set by up to 
~20% in the equivalent width of the H/? absorption line and by 
< 20kms"' in the mean LOS stellar velocity (Fig.[T]). We also 
studied other possible pairs of SSP models, but they were less 
successful in reconstructing the H/3 absorption wings, the metal- 
lic lines, or predicting the Hfi absorption equivalent widths that 
differed significantly from the Tremonti set. Some SSP model 
combinations (such as those of ages 1 1 Gyr and 600 Myr) proved 
useful for certain objects, but the best-fit solution did not differ 
significantly from the final pair of SSPs (11 Gyr and 100 Myr), 
and did not provide optimal solutions for the remainder of the 
sample. 

The best-fit model stellar template for each observed spec- 
trum produced from a linear combination of selected SSPs 
was determined by using the penalised pixel-fitting code of 
Cappellari & Emsellem ( 2004| , and a careful masking of emis- 
sion lines. For the purposes of inferring both stellar kinematics 
and stellar population ages reliably, we increased the signal-to- 
noise ratio by rebinning the FOV with equi-mass Voronoi tessel- 
lation ('Cappe llari & Copin|2003| into irregular compact cells of 
uniform S/N - 50 across the FOV. For other applications em- 
ploying the stellar model to correct the gas emission, we mod- 
elled all the individual observed spectra, without losing spatial 
resolution. 

3.2. Emission lines 

The two observed spectral domains MRl and MR2 contain 
two B aimer lines (Ho' and H/3) and five forbidden-line dou- 
blets of metal ions that represent some of the most fun- 
damental spectral diagnostics for deriving physical parame- 
ters of the ionised gas: [Oiii] ^^4959,5007, [Ni] 198, 5299, 
[Oi] /1/16300, 6364, [Nn] ^^6548,6583, and [Sii]/l/16717,6731. 
In addition to the intensity in individual lines, we derive their 
mean line-of-sight velocities and line-of-sight velocity disper- 
sions and plot their two-dimensional distributions. Not all of the 
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Fig. 1. NGC 2273, comparison of SSP fits in the central bin 
and a bin close to the FOV edge - a young stellar population 
(100 Myr, cyan), an old SSP (11 Gyr, green), an optimal fit re- 
sulting from the combination of two SSPs: 1 1 Gyr and 100 Myr 
(red), and an optimal template based on the Tremonti set of SSPs 
(blue). Spectra binned by Voronoi tessellation at S/N = 50. 



targets were observed in both spectral domains (see Table[2|. 
Where the data are available, we computed the electron density 
from the [S n] doublet and estimated the interstellar reddening 
from the ratio of B aimer lines ( |Osterbrock| 1 989| ) . 

We modelled all narrow emission lines using single Gaussian 
profiles, even though asymmetric lines were expected at many 
locations. Cases for which we obtained large FWHM values 
of Gaussian models led us to identification of spatial regions 
where multi-component fitting might be necessary. We used the 



fit/spec line-fitting software (Rousset 1992| , which has a core 
in C and interface in ESO MIDAS, and permits automated fitting 
over entire datacubes, with constraints applied to the line inten- 
sities, intensity ratios, velocities, and line widths. We restricted 
the number of free parameters by fixing the ratios of the line in- 
tensities within the doublets to comply with the rules of atomic 
physics: 

0.5 < [N I] ^5198/ [N I] ^5200 < 1 .5; 
[O III] ^5007 / [O III] /14959 = 2.88; 
[Oi] /16300/[Oi] /16364 = 3.0; 
[N II] ^6583 / [N II] .16548 = 2.96; 
0.35 < [S II] ^6717 / [S II] .16731 < 1.5 . 

To constrain line parameters more effectively in the auto- 
mated fitting process, especially for the weaker or blended lines, 
the emission lines were associated to form systems of assumed 
common kinematics. The mean LOS velocity and LOS velocity 
dispersion were thus determined for each of the systems, rather 
than the individual lines. The definition of our kinematic systems 
was physically and experimentally motivated. The kinematical 
properties of O^^ ions are known to differ from those of the 
low-ionisation species, due to the possibly different processes re- 
sponsible for the ion formation. The first ionisation potential of 
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elements such as oxygen, nitrogen, and sulphur is almost iden- 
tical to that of hydrogen, and the ions can therefore be formed 
in the same conditions. And indeed, the observations show that 
kinematic properties of the neutral and low-ionisation species 
are usually similar to those of hydrogen. With this justification, 
we associated the weak and blended [N i] doublet with Hfi in 
the same system, while [Oiii] was treated independently. All of 
the MR2 lines (Ha, [Oi], [Nii], and [Sii]) were considered in 
a common system (where possible) to assure a good quality fit 
to lines that are rather weak (mainly [Ni], [Oi], and [Sii]) or 
blended ([S ii] doublet lines, [Nii] with He). 

We corrected the measured FWHM for the instrumental 
dispersion, which was estimated by fitting telluric lines to be 
(7 + 0.5) A. A Gaussian approximation to the line profiles was 
used, hence with a quadratic subtraction of the measured total 
and instrumental dispersions. The resulting dispersion contains 
an error of 0.6 A due to wavelength calibration. The differences 
encountered between the resultant velocity dispersion measured 
in Ha and in H/3 can be almost entirely attributed to their differ- 
ent data quality, far higher for He Where available. Ha results 
should therefore be considered more reliable. 



3.3. Error computation 

Uncertainties in the measurements of the mean LOS velocity, 
which are a function of the signal-to-noise ratio, were studied 

42004), who fitted re- 



for OASIS data in detail by Ferruit et al 



peatedly artificial emission-line spectra with added noise in 500 
Monte Carlo realisations for each S/N value (see also |McDermid| 
et al.|2006 1. Their results were applied in the way recommended 
by the authors in the present paper We show the corresponding 
error bars in one-dimensional plots of mean LOS velocities in 
Figs.[3]-|25] The effects of fitting single Gaussian functions, and 
of observed line-profile asymmetries were not included in this 
treatment. Equally, errors introduced by the subtraction of the 
stellar component, especially in the hydrogen lines, were not in- 
cluded in the error bars. An error of 16 km s ' due to calibration 
was added to the computed error in the mean LOS velocity. 

The random errors in line flux determination were estimated 
with the assumption of Gaussian line profiles, and were com- 
puted as the product of noise in the vicinity of the lines, and the 
line FWHM. This procedure does not include possible system- 
atic errors due to the stellar template determination, which are 
difficult to quantify, as discussed in Appendix B3 of [Emsellem 



et al. ( 2004 ). We plot error bars in the scatter plots of the derived 



quantities such as the line ratios and the one-dimensional cuts of 
electron densities, for all of which Gaussian error propagation 
was applied. For electron densities, the errors were computed 
assuming further that the [S ii] line shifts and FWHMs corre- 
sponded to those of Ha, i.e., we did not consider the effects of 
line blending. 



4. Results for individual objects 

4.1. General comments 

We present colour-coded maps of the line emission intensity, the 
mean line-of-sight velocity, and the LOS velocity dispersion in 
the detected emission lines for each of the Seyfert 2 galaxies in 
our sample. Thanks to our stellar-population modelling, the gas 
velocity maps can be compared with those of stars in the objects 
where the stellar data quality was sufficient for reliable mod- 
elling of the absorption-line profiles. The stellar modelling also 



yielded maps of the mass fractions of the young and old stellar 
populations. 

Further information about the objects can be inferred from 
the ratios of the emission-line intensities. We derive electron 



density from the [S ii] line ratio, as described in [Osterbrock 
( |I989| l, by assuming a temperature of 10 000 K (our data do not 
contain a direct temperature diagnostic). We estimate the dust 
distribution from the Balmer line ratio. The Ha/H/3 should have 
a fixed intrinsic value depending on the physical conditions in 
the environment where the recombination lines are formed. We 
assumed the intrinsic value of 3.1, which is typical of NLRs 
( [Oste rbrock 1989). By assuming any excess of the measured ra- 
tios to be due to dust along the line of sight, we computed the 
interstellar reddening E(B - V). 

The ratios of forbidden-line to recombination-line intensities 
are a measure of the ionisation and excitation state of the gas, of 
its density, temperature and metallicity, as well as the nature of 
the ionizing source. We plot maps for each of the standard diag- 
nostic ratios [O m]/l5007/Hy6 (excitation map), [Oi]/16300/Hq', 
[Nn]/i6583/Ha, and ([S ii]^6717-h[S ii]^6731)/Hq'. More infor- 
mation is conveyed by the spectral-diagnostic diagrams defined 



by |Baldwinet al. (1981); and Veilleux & Osterbrock ( 1987), but 
in our case, these diagrams are plotted for spatially resolved data 
for each observed galaxy. The diagrams provide further possi- 
bilities for future detailed modelling of the ionisation structure, 
similar to that performed e.g., by Bennert et al. (2006b) , or Allen] 
|et al. ( 1999| l for long-slit spectroscopic data. 

In our diagnostic diagrams, we compare the different existing 
classification schemes for emission objects commonly used in 



the literature, including the recent development (Kewley et al. 
2001 ; Kauffmami et al.|2003[|Stasiriska etal. 2006. Kewley et aT 



-^i I 1 1 I 1 1 : I 

2006|l, as well as the traditional schemes, some of which were 



based on a single pair of lines and were reviewed by Ho et al 
( |1997) l. However, we note that the models correspond to global 
spectra, and one must be careful when interpreting the spatially 
resolved diagrams (see ,Stasiriska,2009) , which require detailed 
ionisation modelling for a proper understanding. Therefore, the 
plotted curves that correspond to the afore mentioned models 
play an informative role only. 

All the plots presented in Figs.[2)-[25l show the results of 
single-Gaussian modelling of the emission lines. This simpli- 
fication is clearly insufficient for a precise reconstruction of 
many spectra, where two or more components are present in 
each emission line, as specified below for the individual objects. 
The single-component fitting obviously limits the interpretation 
of the presented data, although the results provide a guide to 
the location of important structural and kinematic features in 
the objects. The applied modelling approach simulates the ef- 
fects of lower resolution spectral data with unresolved multi- 
ple components. The resulting intensity in the lines corresponds 
to the integrated intensity of all components or to the inten- 
sity in the strongest component, depending on the line shape 
and the fit. Similarly, the fitted line-of-sight velocities trace ei- 
ther the motion of the strongest component at a given loca- 
tion or a "weighted average" of the components. The velocity- 
dispersion maps serve as efficient indicators of locations of the 
multi-component spectra. While it would be impossible to in- 
spect all spectra individually due to the enormous amounts of 
data, the large-dispersion areas of the FWHM maps trace the 
regions of complex emission-line profiles, often in unexpected 
locations, not necessarily at the galactic centres. 

We display profiles of the mean LOS velocity and electron 
density along the major axes of emission, and scatter plots of ra- 
dial variations in the diagnostic line ratios as complementary in- 
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formation. Most of the data presented were clipped at S/N = 5, 
with the exception of the weak lines [Ni] and [Oi], for which 
only at a few spatial pixels in the centre the signal-to-noise ratio 
exceeded the value of two, and therefore we applied clipping at 
S/N - 1 to plot the two-dimensional maps. The surface bright- 
ness is in the units of 10"'^ J s"' arcsec"^ m"^, the mean LOS ve- 
locities, and LOS velocity dispersions in kms '. The velocities 
were corrected for heliocentric motion using the rvcor routine 
of IRAF, and plotted relative to the systemic velocity of each 
galaxy, which was taken from NED (see Table[T]). 



We compared the measured fluxes of emission lines with 
those determined by |Falcke etal.| ( 1998 1. The authors provide the 
fluxes in [Oni] and HQrH-[Nii] obtained from HST observations 
of Mrk 34, Mrk 348, and NGC 4388, using circular apertures of 
two different sizes for each object. We simulated the apertures by 
summing the spectra over circular regions centred on [0",0"] 
and of the same radii as in the cited paper The emission-line 
fluxes were measured for spectra that had a stellar component 
modelled by a combination of two SSPs (11 Gyr and 100 Myr) 
from [Bruzual & Chariot] ( |2003[ ), as described in Sect. |3.1| and 
with no emission-line model applied (to avoid introducing fur- 
ther sources of error). Results of the comparison are summarised 
in Tabled 



Agreement between the two data sets is closer for larger 
apertures, as expected. The flux that we measure in regions of 
radii smaller than 1" is affected too much by the spatial PSF 
(dominated by ~ 1" seeing), which is an order of magnitude 
higher than that in the HST observations and redistributes much 
of the observed flux outside the selected area. Most of the fluxes 
that we obtained for apertures of sizes > 1" are 10% - 25% 
lower than previously published values. Possible reasons for the 
discrepancy include flux-calibration errors, errors acquired dur- 
ing data reduction or due to an imprecise determination of the 
stellar-continuum level. However, it is impossible to draw con- 
clusions about systematic differences based on only three cases. 



Table 3. Comparison of our emission fluxes with those of Falcke 
|et al.| ( [T998[ l. OASIS values correspond to fluxes corrected for absorp- 
tion usmg a combination of SSPs as described in Sect. 3.1. Units are 
lO^^Wm^l 



Galaxy 


Aperture 


OASIS data 


'Falcke et al.'(1998^ 




radius 


[Oiii] 


Ha+[Nii] 


lUiiiJ 


H(H-[N n] 


Mrk 34 


2" 


4 800 


3 400 


6100 


3 600 




0.5" 


500 


300 


1300 


700 


Mrk 348 


1.5" 


3 000 


2 500 


4100 


2 000 




0.4" 


1500 


1000 


3 200 


1500 


NGC 4388 


6" 


9 800 




1 1 700 






I" 


3 200 




3 900 





The following subsections are devoted to the description of 
the individual objects, the observed features, set into the larger 
multi-wavelength context found in the literature. We note that 
we will not be restating the wavelengths of the emission lines 
that were listed in Sect. 3.2 except the cases where confusion 
would be possible. 



4.2. Individual objects 
Mrk 34 

Mrk 34 is known to be a luminous FIR source with a power- 
ful water maser ( Henkel et al.|2005 l. Optical line emission was 
found to be spatially correlated with radio continuum emission, 
PA ~ 150° - 165° (|Nagar & Wilson |[T999HFalcke et al.||1998l 
lUlvestad & Wilson] 1984a| ): 

Mrk 34 is the most distant object in our OASIS sample 
and in this case our FOV corresponds to the central ~ 5 kpc. 
Because of the high redshift of the galaxy, the [S ii] doublet is 
located outside the observed spectral range, and we are there- 
fore unable to infer all of the spectral-diagnostic ratios, and the 
electron density. The observed emission in all of the observed 
lines has a similar morphology, extended along the direction of 
PA ~ 140° - 150° (Fig.|2]). The excitation map [Oiii]/Hj6 has the 
same major axis as the [Oiii] emission, and more of a cone-like 
shape. The other diagnostic line ratios, however, have low values 
in the direction of the cones, while the maximum values are dis- 
tributed along the minor axis (Fig.|2]i. Reddening computed from 
the HafH/3 ratio possesses a maximum in the southern part of the 
FOV, ~2 kpc from the central maximum. The spatially resolved 
diagnostic diagrams indicate that the emission corresponds to a 
strong Seyfert regime. 

The map of the mean LOS velocities in Ho- (Fig.|3]l provides 
an example of twisted S-shape isocontours. They are not evi- 
dent in [Oiii] and H/3, which may partially be caused by the 
smaller FOV in MRl. No galactic bar or other signature of non- 
axisymmetric potential has been reported so far, and the S-shape 
might be the first piece of kinematical evidence of departures 
from axial symmetry. The low quality of the stellar data did not 
allow comparison of the stellar and gas velocities. The object 
is characterised by multi-component emission-line profiles in a 
region of elongated shape perpendicular to the emission major 
axis, as visible in the velocity dispersion map. 

Mrk 622 

The radio emission in this galaxy at the 20 cm wavelength is 
extended along PA - 0° ( Nagar et al.|1999 1. No prominent axis 
in either the optical emission or the excitation map was detected 
on the scale of 5" by Mulchaey et al. ( 1996a), while Schmitt 



et al. (2003a I reported elongation in [Oiii] emission at PA = 
55° in the inner 1", perpendicular to the host-galaxy major axis. 
|Shuder & Osterbrock| ( 1981 1 discovered multi-component [Oiii] 
profiles of total FWHM - (1050 + 50)kms in contrast to 
the nearly-Gaussian shapes of the Ha, [Nii], and [On] lines of 
FWHM = 350 ± 75kms They concluded that two emission- 
line regions of different ionisation are present, with estimated 
temperatures of ~ 1 1 000 K in [O iii] and ~ 6000 K in Ha, [N ii] , 
and [On]. 

The galaxy is among the most distant in our sample, which 
has an effect on the spatial resolution that we attain. The maps 
of line emission intensities do not show any particular direction- 
ality, except for the Ha and [N ii] emission, which has hints of 
extending both to the east and the south (Fig.[ 4ll, in ag reement 
with the optical and radio observations of [Schmitt et al. ( 2003a| l; 
and Nagar et al. ( 1999[ l. The [Oiii] lines are highly asymmet- 
ric and their multiple components are detected throughout the 
field of view (Fig.|6]), with a separation between the components 
as large as 600 kms"', not present in the other lines, confirm- 
ing the results of |Shuder & Osterbrock| ( [T981| l. Single-Gaussian 
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modelling of [Oiii] thus results in large-FWHM fits, as seen in 
Fig.|5] 

The one-component fits have a strong impact on the interpre- 
tation of the velocity field of this object. The relative intensity of 
the two line components in [O m] varies across the FOV, with the 
blue component becoming stronger in the east and the red com- 
ponent in the west. The fitted one-component Gaussian model 
of the line profile is always biased towards the stronger compo- 
nent. As a result, our velocity map in [O in] is not a good ap- 
proximation of the true motion and a detailed multi-component 
modelling of line profiles is necessary. We assume that this is 
the major reason for the dissimilarities between the hydrogen 
and [Oiii] kinematics presented in Fig.|5] While the apparent 
kinematic axis of [Oiii] lies along PA ~ 90°, that of Ha (to- 
gether with the low-ionisation lines) is perpendicular to it, at 
PA ~ 180°, measured in the outer parts of the Ha map, the cen- 
tral regions being characterised by S-shape isovelocity contours. 
The stellar velocity field, as far as can be deduced from our data, 
has a kinematic axis at PA ~ 120°, and it is impossible to de- 
cide whether the twisted contours are present due to the low data 
quality. 

The data are of low quality and their derived velocities do 
not agree with those measured from Ha. We note that velocity 
colour scales in Fig. [5] were not homogenised to the same lim- 
its, due to the problematic fitting results. Two-component line 
profile modelling is necessary for [O ni] lines, as the ambiguous 
interpretation of present results demonstrates. 

The electron density, the interstellar reddening, and the gas 
excitation all reach their maximum values at the galactic centre. 
In contrast, the diagnostic ratios of line intensities, except for 
[O iii]/HyS, reach their minima in the central region. Interestingly, 
the spatially resolved diagnostic diagrams indicate low excita- 
tion across the entire FOV (Fig. [4]). The observed spectra pos- 
sess strong Hy6 absorption wings, which affect both the selection 
of the stellar population model and the predicted absorption in 
HyS, which may be overestimated. 

Mrk1066 



WTiile being classified as a marginal Seyfert by Osterbrock & 



Dahari ( 1983 a high-excitation component of the spectral Imes 



was found by 
at PA ~ 135 



Bower et al. ( 1995 1. Radio emission was observed 
by'Haniff^ et al.|(|1988); and|Ulvestad & Wilson 



citation, and a small velocity dispersion (Figs.|7]and[8]l. We con- 
clude that it is probably a region that has not been ionised by the 
AGN. Excitation is highest in the central parts and extends asym- 
metrically to the north-west side along the major axis of [O iii] 
emission, as well as along the minor axis to both sides from 
the centre. Interestingly, the interstellar reddening extends asym- 
metrically along the major axis in the opposite direction, to the 
south-east (Fig.|7]i. The maps of diagnostic ratios [Oi]/Ha and 
[N ii]/Ha exhibit surprising distributions of high values within 
an extended region along the minor axis (Fig.]?}. The diagnostic 
diagrams are consistent with the interpretation of the object as a 
marginal Seyfert (Fig.[7]i. With a few exceptions, all of the spec- 
tra are characterised by an [O iii]/HyS intensity ratio that is below 
the empirical Seyfert limit, including the central parts. However, 
a proper disentangling of the asymmetric lines may confirm the 
high-excitation component ((Bower et al. 1995)). The electron 
density reaches peak in the central region, the maximum being 
off-centred by -200 pc (Figs.[8]and|9]l. 

The Ha velocity map (Fig.[8]l displays a clear S-shaped pat- 
tern. Velocities in [O iii] exhibit important differences from the 
hydrogen data, seen mainly in the one-dimensional cuts. We 
found no multi-component line profiles for this object (Fig.|9]l, 
even though we observe a large-FWHM region extending along 
the minor axis, which is usually associated with line splits in 
other objects of the sample. On the other hand, most lines have 
blue wings across the entire FOV. 



NGC262 (Mrk 348) 

The host-galaxy is orientated almost face-on and has a close 
companion, NGC 266. The nuclear region has been a popular 
target of observations across the entire electromagnetic spec- 
trum, because of the difficulties to interpret orientations of in- 
dividual components observed in the different spectral regions 
within a simple model ( [Anton et al. 12002) 1 . Capetti et al. ( 1996^ 
reported a prominent central dust lane of scale height 50 pc, 
perpendicular to the optical [Oiii] emission in an HST image. 
Evidence of a torus or related structure was found in IR im- 
ages by Simpson et al. ( 1996). A polarised broad Ha component 



1989| l, and the narrow-band images indicated an alignment of 



(Miller & Goodrich 1990[ l and a hard X-ray emission (Awaki 



|et al. 1991) inferred a hidden Seyfert 1 nucleus. The radio emis 



the optical [Om] emission with the radio jet ( Bower et al.|1995 



Mulchaey et al.| 1996a i. The direction of both types of emission 



sion is in the form of a triple radio source of spatial extent 0.2" 
along PA = 168° ( |Neff & de Bruyn||T983] l, variable on the 
scale of months. Megamaser emission detected in the central few 
then agrees approxi mately with the orientation of the galactic milha rcseconds (|Falcke et al. 2000; Xanthopoulos & Richards 
bar of PA ~ 143° (Mazzai-ella & Boroson 1993; Bow er et al] '2001; Pe ck et al.| |2003) was interpreted as an interaction be- 



|1995) . According to Bower et al. (19 95 ), the velocities in Ha 
[N n] are consistent with rotation, probably in the galactic disc, 
while [O m] shows little or no rotation, which they interpret as 
an outflow. No signs of hidden BLR were inferred from spec- 
tropolarimetry ( Miller & Goodrich| |1990). No central dust lane 
in high-resolution HST images, nor other evidence of a torus 



(such as ex c ess reddening at the c entre) was found by Bower 
et al.| ( |1995] ). [Zhang & Wang| ( |2006| l attributed the non-detection 
of a hidden BLR to orientation effects and/or the low mass of the 
central black hole. 

The emission that we observe in all the lines is slightly ex- 
tended along PA ~ 130° - 140°, i.e., parallel to the galactic bar 
A prominent emission feature located in the east is seen in the 
hydrogen lines and partially in [Nii] and [Sii], while absent in 
[O III], and is associated with high interstellar reddening, low ex- 



tween the radio jet (situated close to the plane of the sky) and 
a molecular cloud. Optical [O iii] emission is well ahgned with 
the radio on sub-arcsecond scales ( [Falcke et al.[[1998) [Schmitt[ 
et al. 200 3a|, while on the arcse cond scales it is extended along 
PA ~ 185° ( (Schmitt et al.[2003al l, with a secondary blob of emis- 
sion ~ 0.9" south of the nucleus. A hint of ionisation cones was 
found by [Simpson et al.| ( [1996[ ); and [Mulchaey et al.[ ( [1996a| l 
in the ionisation map, but the evidence is marginal according to 



Anton et al. (2002 ). NIR emission map obtained with an IFU by 



Sosa-Brito et al. (2001) does not have any prominent direction 
on the scale of 3", although contours in the inner parts seem to 
have a slight directionality along PA ~ 10° and PA ~ 90°. A ring 
of H n emission on the scale of ~ 1 kpc was reported by "Anton] 



let al.((|2002); and Mulchaey et al. (1996a). The galaxy possesses 



a giant H i envelope on the scales of hundreds of kpc ( Heckman 
ietal.|1982| ). 
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The face-on orientation of the galaxy is confirmed in our data 

the minimum and 



by the low stellar LOS velocities (Fig. 1 1 



maximum values that we obtained from our fits are separated 
by 70kms"' only. We estimate the kinematic stellar axis to be 
at PA ~ 70°. This Seyfert galaxy has one of the lowest contri- 
butions from young stars of the entire sample. Mass fractions 
of young stars above several percent have been identified in the 
central ~250 pc only. 

The emission maps of most of the emission lines are slightly 
extended along PA ~ 10° (Fig. 10 1, which is consistent with 
Sosa-Brito et al. ( 2001| l, and differs by approximately 20° from 
the radio jet reported in the literature. We partially resolve the 
secondary emission maximum known from the literature, ~ 1" 
south-west of the centre. The gas excitation is high, and all the 
regions are high in the Seyfert regime according to the spa- 
tially resolved diagnostic diagram (Fig. 10 1. The regions of high- 
est excitation are situated along the major axis of emission, 
within a PA range of ~ 5° - 30° . The [Oiii]/H/3 ratio decreases 
steeply with radial distance, as seen in the corresponding graph 
of Fig. 1 1 We note that the ratios of low-ionisation lines to Ha 
show prominent minima south-west of the galactic centre, in the 
region of the receding velocity maximum (see paragraph below). 
In addition, [S ii]/Hci' has a minimum at the centre, surrounded 
by a patchy ring of high [S ii]/Ha on scales of ~600pc, which 
corresponds to a ring of low electron densities derived from the 



on the scale of ~30pc ( |Kukula et al.||1999) l. Stellar-population 
synthesis modelling of long-slit data of 3 A resolution performed 
by R aimann et al.| ( )2003) revealed mostly intermediate-age stel- 
lar populations (~ 1 Gyr) in the central kpc and an increasingly 
dominant old population (~10Gyr) outside this radius. 

Our OASIS data for this object were obtained as a mosaic 
of three partially overlapping fields. We confirm strongly colli- 
mated emission in [Oiii], Hii, and [Nii] extending to radii of at 
least ~3 kpc (Fig.[T3]l. In the other observed lines, the S/N is in- 
sufficient to obtain reliable data to these distances. We found a 
misalignment between the emission in [O m] (PA ~ 65°) and the 
low-ionisation species (PA ~ 85°). The velocity map in Ha sug- 
gests circular rotation. As in most of the targets of the OASIS 
sample, we found a region of large FWHM, stretching along 
the minor axis (Fig. [14]). The region is characterised by multi- 
component line profiles, usually in the form of red or blue wings 



[Sii] ratio (Figs. 10 and 11 



in [O III] and split line profiles in low-ionisation species (Fig. 15 i. 
Outside the central ~600 pc, the emission lines are narrow. 

Our stellar data were not of a high quality for this object. The 
maximum contribution of young stars was fitted in the nuclear 
region. The map of stellar velocities could not be interpreted re- 
liably. 



NGC 2273 (Mrk 620) 



The gas velocities (Fig.[TT) are the most surprising detec- NGC 2273 is a barred galaxy ([van Driel & Butel[T99T), with 



tion. We observe two spots of velocity maxima of opposite sign, 
aligned at PA ~ 25°, i.e., approximately parallel to the major 
axis of optical emission. The two spots of > 100 km s"', present 
in [Om] as well as in low-ionisation lines, about 200 pc in size, 
are placed symmetrically ~300 pc from the emission maximum. 
We assume that the observation corresponds to a rotating ring, 
inclined with respect to the galactic disc. However, there is a 
clear asymmetry in the radial profile of the two velocity peaks 
(Fig. 



star-forming nuclear ring ( Permit et al. 2000[ ). Deviations 



from axial symmetry were confirmed by GMOS IFU observa 
tion of the stellar kinematics in the central regions (Barbosa 
2006[ l. Radio emission extends along PA ~ 20° at the 
20 cm wavelength, and two components along PA ~ 90° are 
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The maxima of the LOS velocity dispersion in Ha are lo- 
cated at the positions of the velocity peaks, which is not the case 
for [Oiii]. The low-ionisation emission lines also have signifi- 
cantly asymmetric profiles across the FOV, unlike [O iii] where 
mostly blue wings are present (Fig. 12 1. The low-ionisation lines 
are characterised by prominent blueshifted and redshifted com- 
ponents shifted by 300 - 400 km s"' from the main component. 
The red asymmetries were found mainly in the southern part of 
the FOV. The vicinity of the south-western spot of maximal ve- 
locities is characterised by rapid spatial variations in the line 
shapes and switching between red and blue additional compo- 
nents. A careful disentangling of the line profiles is necessary 
before further conclusions can be made. 

We detect another systematic motion in [Oiii], at the edge of 
the FOV, of radii ~1 kpc, where we observe velocities with a po- 
larity that is reversed with respect to the nuclear ring. We verified 
that it is not the effect of noise and suggest that it might corre- 
spond to the ring identified by Anton et al., ( ,2002) ; and |Mulchaey, 
|etaL| ( |1996a| l. 

NGC 449 (Mrk 1) 

The galaxy is well known to possess a prominent water maser 
in its nucleus ( |Braatz et al.||1994) l. No evidence of a hidden 
Seyfert 1 nucleus was found, in either polarised light (Kay' 1994) 
or in the infrared ( Veilleux et al.|1997] ). Mulchaey et al. (1996a) 
reported [O iii] emission extending to ~ 3 kpc, co-aligned with 
both the galaxy major axis, PA - 83° and 20 cm radio emission 



evident at the 6 cm wavelength in the central 1" (lUlvestad & 
I Wilson|1984bl|Ho & Ulvestad|2001| l. HST observations resolved 
a jet-like [Om] structure extending 2" east of the nucleus, at 
PA ~ 90°, aligned with the double radio source (Ferruit et al. 
20001. Mulchaey et al. ( 1996a[ l reported extended [Om] emis 



Moiseev et d?] ( |2004| i 



sion reaching to -6.5" at PA = 134° ± 10° 

identified important differences between the velocities in [O m] 
and the Balmer lines. 

Our data for this object are limited to a single spectral do- 
main (MRl). While the [Om] emission does not exhibit any 
prominent axis (only a slight elongation is visible along PA ~ 
95°), Hfi emission seems to follow the direction of PA ~ 20°. 
These differences in behaviour are emphasised in the [O iii]/HjS 
map, where extended cone-like features appear along PA ~ 100° 
(Fig. [16]). The velocity field in [Om] lines shows S-shaped iso- 
contours (Fig.[T6|, which might indicate deviations of the grav- 
itational potential from axisymmetry. The data quality in Hfi is 
insufficient to confirm or exclude a similar velocity pattern in 
hydrogen. However, one-dimensional cuts of the velocity fields 
show important systematic differences between LOS velocities 
in [O III] and H/3. The stellar velocity field is well measured for 
this galaxy and shows rotation in a direction similar to gas, with 
the kinematic axis at the outer parts of the FOV approximately 
parallel to the H/3 photometric major axis. No asymmetric or 
multi-component emission-line profiles have been found, apart 



from weak, blue wings (Fig. 16 1. 



NGC 2992 

The host galaxy disc is highly inclined at / ~ 70°, which makes it 
suitable for extraplanar emission studies. The galaxy has a com- 
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panion NGC 2993, which is located ~3' south-east, connected 
with a tidal bridge, and whose tidal forces might have an im- 
portant effect on the kinematics. Radio observations at 20 cm 
found a 25" structure along the major axis of the galaxy, and 
a 90" one-sided extension at PA ~ 100° - 130°, close to the 
galaxy minor axis ( [Ward et aL 1980| l. The 6 cm radio emission is 
concentrated in the central parts, forming a "figure-eight" struc- 
ture at PA = 160° (Ulvestad & Wilson 1984b; Wehrle & Morrisl 



1988| l. Soft X-rays are detected up to scales of 35" - 45" along 
the galaxy's minor axis, and are co-spatial with the 20 cm radio 
emission. The [O iii] emission forms two sharp-cut cones with a 
projected opening angle of 120° (Marquez et al. 1998 ), their axes 
aligned with the galactic minor axis (PA - 120°) and approxi- 
mately with the radio structure. A disturbed dust lane was iden- 
tified at PA ^ 1 5° byJMarquez et al.| ( [T998 ) and at PA = 30° by 
Veilleux et al. ( |2001 1. Multi-component optical emission lines 
were detected, and the velocity field was suggested to be a super- 
position of rotation and an outflow (Marquez et al.,1998| Allen 
let al.|1999l[Veilleux et al.\2Wl\ . 

The FOV in OASIS observations enables only one of the 
emission cones to be imaged, located south-east of the dust lane. 
The emission structure that we observe (Fig. 17 1 is elongated at 
PA ~ 30° in all of the available emission lines, parallel to the 
dust lane observed by Veilleux et al. ( 2001) ). The emission bends 
in the southern part towards PA ~ 210°, especially in [Om], and 
to PA ~ 45° in the north of the FOV. [Om] has a significant 
emission extension in the eastern part of the FOV, associated 
with high receding velocities, while the low-ionisation lines ex- 
hibit a blob of emission to the south-west, associated with high 
extinction. 

The maximum of the emission in all lines coincides with the 
maxima in the electron density, the interstellar reddening, and 
the velocity dispersion (with the exception of [O m] where the 
FWHM maximum is shifted by ~l" to the east, see Fig. 18 i. At 
the same time, the maximum emission region is characterised by 
local minima of diagnostic line ratios (Fig.[T7|). 

The LOS velocity fields (Fig. 18 1 obtained from different 
emission lines are similar, characterised by receding motions 
in the northern part and approaching in the south and east of 
the FOV. Comparison with the stellar velocities (with their kine- 
matic axis in the NE-SW direction) suggests that the gas has a 
rotational component aligned approximately with stars. Despite 
the apparent kinematic similarities, the MR2 spectra had to be 
considered to correspond to two kinematic systems, by fitting 
the [Sii] and [Oi] kinematics independently of Ha and [Nii], 
while the unified fit was impossible. The grouping of the emis- 
sion lines into the two systems was a compromise between the 
well-determined positions of lines and the number of parame- 
ters needed to place enough constraints on the weak lines and 
the blended lines. We illustrate the blending of the Ha and [N ii] 
lines in the spectra of Fig. 19 Outside the central ~2" (350pc), 



all emission lines are however naiTow and non-blended. The dif- 
ferences in the motion of the individual ionic species are de- 
picted in the one-dimensional cut along the major axis (Fig. 18 1: 
the velocity curve in [O i] and [S ii] is shallower (by ~100 kmT"^) 
than in the other lines, whereas the most prominent velocities are 
observed in [Om], especially in the northern part. The velocity 
curves that we obtained by one-dimensional cuts along diff'erent 
axes are consistent with long-slit spectroscopy by Marquez et al. 
( [T998l l. 

We detect double-component profiles of emission lines, par- 
ticularly for [O m] and H/3, while the other lines are too affected 
by either blending or noise. We show the asymmetric MRl pro- 



files in Fig. [19] in a combined plot of spectra from three loca- 
tions. Broad wings of Balmer emission lines were detected in the 
central ~2" (~350pc). However, their relation to BLR emission 
(already suggested e.g., by Ward et al. 1980; B oisson & Durret| 
1986) can only be confirmed after a careful disentangling of all 
the line components. 

Our stellar population modelling revealed spatial variations 
in the contribution of young stars, which is highest in the 
strongly emitting region close to the dust lane, while mostly old 
populations are present in the rest of the FOV (eastern part). 

NGC 3081 

The host galaxy is an early-type spiral with a weak large-scale 
bar and a nuclear bar (! Buta|[T990 l), and four well-defined reso- 
nance rings: nuclear, inner and two outer ones (Buta & Purcell 



1998). The brightest emission-line region inside 2" has a Unear 
morphology, well aligned with the radio emission ( Nagar et al. 



1999). Slitless spectroscopy by Ruiz et aL| ( 2005 ) revealed two 
major kinematic components in the [Om] line at — 250 km s ' 
and -i-50kms"'. 

We observed the object in one spectral domain (MRl) only. 
The photometric axis in both [O m] and H/3 stretches along the 
PA ~ 155° - 165°, even though it is less pronounced in [Om] 
(Fig. [20]). Unlike in [O m], the emission in H/3 is elongated asym- 
metrically to the south-east. 



The mean LOS velocity fields of [Om] and H/3 (Fig. 20 1 



show similar structures, forming an S-shape in [Om]. However, 
important differences between [O iii] and Hy6 velocities are en- 
countered in one-dimensional cuts of the velocity maps. The S- 
shape is not reproduced in the stellar velocities, derived from 
high-quality absorption-line data in this object, which have a ro- 
tational pattern with a kinematic axis PA ~ 90°. The anomalous 
values of stellar velocities detected in several bins of the map in 
Fig.|20] do not correspond to true velocities, and are results of 
poor fits in these particular positions. Interestingly, the northern 
region of receding [O m] velocities is characterised by minimum 
values of the line width. We found no prominent asymmetries in 
emission line profiles (Fig. 20 1. 



Our stellar population modelling revealed low contributions 
from young stars, compared to the other Seyfert nuclei in the 
studied sample (Fig. 20 ). Nevertheless, the SSP ages form a clear 
spatial structure with a boundary at ~ 800 pc separating older 
populations inside this radius and younger outside. 

NGC 4388 

The galaxy belongs to the Virgo cluster ( [Phillips & Malin|1982t 



Yasuda et al. 



78° 



T997|). The galact ic disc of PA ~ 90° is inclined at 
leux et aL]|1999a| ), with the north rim on the near 



side. The nucleus is a strong X-ray source (e.g ., [Hanson et aL 
1990; Ta kano & Koyama|[T99T] [Lebrun et"aL]|1992[ [Iwasawa 
et al.,1997[ ). Radio-emission morphology suggests a collimated 
AGN-driven outflow reaching ~ 200 pc south of the nucleus 
(e.g., |Stone et al.|1988[[Falcke et al.TT99 8 ). Weak broa d Hg emis- 
sion impUes that a hidden type-1 nucleus is present ( [Filippenko 
|& Sargent||1985[ ). Extended optical emission has two compo- 
nents of different excitation, one associated with the galactic 
disc and another reaching 50" (4 kpc) above the galactic plane 
in the form of two opposite cones (e.g., | Pogge 1988 ; ' Corbin| 
let al.||1988l I Falcke et al.|[T998l [Veilleux et al. 1999b ). While 



Rubin et a^ 



(1997 1 argued that anomalous kinematics in the 



central parts of the galactic disc were indicative of presence of 
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a discrete, rapidly rotating circumnuclear disc, | Veilleux et al.| 
( 1999a|b ) found that the gas kinematics were consistent with el- 
liptical streaming in the bar potential, plus a bipolar outflow ex- 
tending out of the galactic plane (and excluded other models for 
the extragalactic component). They argued that the bar-induced 
non-circular motions may account for the emission-line splitting 
of the maximum amplitude of ~150kms"' detected symmetri- 
cally out to radii of ~10" (~1 kpc) along the disc, i.e., close to 
the end of the bar. 



In the case of this nearby galaxy, our FOV is restrained and 
covers the central ~400pc only, which is insufficient for inter- 
preting the large-scale structures known from e.g., Veilleux et al. 
( | 1999b) , but we sample the central region with the spatial scale 
of ~ 80 pc, obtaining thus a highly detailed picture. The maxi- 
mum of emission within our FOV is located in the south-west 
(Fig. 21 1, and the emission intensity isocontours have strongly 
irregular forms, pointing to the north-east PA ~ 30° - 45°, i.e.. 



a direction consistent with the extragalactic component reported 
by Veilleux et al. ( 1999b I. The H/3 morphology differs slightly 
from [O III] in the norther part of the FOV, where it bends by 
~ 90° at ~200pc. The map of [Om]/Hy6 ratio reveals high ioni- 
sation in the region of maximum [O m] emission extending along 
PA ~ 45° (Fig.|2T). 

Due to good-quality stellar data in this object, our map 
of mean LOS stellar velocities exhibits a clear pattern, which 
is similar to that observed in Ho- emission by [Veilleux et al. 
( |1999b[ ). Because of the large FOV covered by the Fabry-Perot 
interferometer observations, they were able to identify a strongly 
S-shaped form in the isovelocity contours in the galactic disc. 
For the gas velocities that we observe, our limited FOV only al- 
lows us to identify agreement between the general orientations 
of the central velocity features in our data and those of VeiUeux 
|etaL] ( |1999b| |. 

We detect two regions of high velocity dispersion, especially 
in [Oiii] maps: in the south close to the peak emission, and in 
the north of the FOV. These regions are characterised by anoma- 
lous emission-line profiles, the former with weak red wings, the 
latter with a marked additional blueshifted component separated 
by ~ 800 - 900 km s"' from the main emission component both 
in [Oiii] and HyS (which is significantly higher than reported 
by Veilleux et al. 1999b) l. Examples of spectra and their one- 



component fits are presented in Fig. 21 from which it is clear 



that the velocity maps follow the strongest component of the 
lines. The additional redshifted and blueshifted components are 
restricted to the two regions and never outshine the main com- 
ponent. 

The stellar population modelhng shows the largest contribu- 
tion of young stars (~ 10% of mass) close to the point of in- 
version of the stellar velocities. The observed field is otherwise 
dominated by old stars, consistent with the findings of |Storch>| 
Bergmann et al.|(|1990|l. 



NGC 5728 

The host galaxy has a bar along PA = 33°, a nuclear bar at PA ~ 
86° - 90° ( |Wozniak et a l. 1995 ; Emse llem et al.|2001| l, and three 
rings. The nuclear ring has dimensions 9" x 7", and the major 
axis at PA = 20° . The Seyfert core is known as one of the most 
spectacular examples of the biconical emission-line structure in 
optical observations, with the north-west cone at PA = 304° and 
the south-east cone at PA = 118°, and opening angles ~ 55° - 
65° ( |Pogge|1989[ [Wilson et al.|1993| l. The emission maxima in 



both Ha;H-[Nii] and [Oiii] are located ~1" from the cones' apex 
dWilson et allT993] l. 

The kinematics of NGC 5728 has been characterised by con- 
troversial conclusions, and so far none of the proposed models 
(radial outflow, inflow or non-axisymmetric rotation in a barred 
potential) has managed to account fully for all the observational 
results. The inflow model is supported by the observed asym- 
metry in the line-emission distribution and the degree of ioni- 
sation ( [Schommer eraL][T988l [Wilson et al.[ [r993; Riff^e l^eTaL] 
120081): the shorter and less ionised north-west cone can be ex- 
plained by obscuration and interaction with the galactic disc if 
the north-west is the far side (see Figs. 22 and 23 for a better 



understanding). On the other hand, HST observations revealed a 
sharp apex to the north-west cone, whereas the south-east cone 
terminates with a blunt end 0.5" from the core, suggesting ob- 
scuration by foreground material ( Wilson et al.|1993) . Together 
with an excellent alignment with the radio emission, this might 
be considered evidence of an outflow ( [Schommer et al.|l"988) . 

Observed with the use of the Fabry-Perot spectrometer by 
Schommer et al. ( [1988 1, strongly S-shaped Ha isovelocity con- 
tours in the central 10" (2 kpc) suggested that these features orig- 
inated in the bar potential. On the other hand, the bar is in a po- 
sition for which models by [Roberts et al. ( 1979|_do not predict 
velocity field distortions, therefore Schommer et al.[ ( "l988| l drew 



no conclusions about the importance of the bar SPH simulations 
by Perez et al., (2004), which included a bar, failed to reproduce 
the observed velocity field and the authors proposed that either 
an additional nuclear bar or radial flows are present. 

The major axes of the emission structures determined from 
our OASIS data are aligned at PA ~ 125° (as measured from the 
outer contours). The emission is asymmetric in the two lobes, as 
reported before in the literature (see above). The north-west lobe 
shows an elongation at PA ~ -40° in the inner 4" (~ 800 pc), 
consistent with the HST observations of Wils on et al.[ ( [1993] l, 
while the more prominent south-west lobe bends towards PA ~ 
1 15°. We partially resolved the secondary maximum of emission 
located ~2" north-west of the maximum, which was reported 
e.g., by [Wilson et al. ( 1993 1. The hydrogen emission maps show 
a prominent arc perpendicular to the lobes, at a distance of ~ 
5"(1 kpc) north-west of the centre. We interpret this as part of 
the ring reported e.g., by Wilson et al. ( 1993 1. 

The gas velocity maps and their one-dimensional cuts 



(Fig. 23 1 show a strong gradient of LOS velocities along the ma- 
jor axis of emission, in the vicinity of the dynamic centre of the 
galaxy. The velocity maxima of opposite signs (approximately 
±300 km s"') are separated by only 400 pc and are located at the 
primary and secondary maxima of emission intensity. A compar- 
ison with the Fabry-Perot observations (Schommer et al.|1988| l 
provides a different view of the velocity maps across a wider 
field, providing thus a larger context which is inaccessible with 
our limited FOV: a strong S -shape is found, which must be taken 
into account in the interpretation (in a future paper). The hydro- 
gen velocity field of the outer parts of our FOV (outside the cen- 
tral 500 pc) has an orientation corresponding to the stellar mo- 
tions that we measured. We suppose that the interpretation of the 
kinematic structure of NGC 5728 will be aided by rigorous dis- 
entangling of multiple-component emission lines that we found 
to be distributed within an elongated region approximately per- 
pendicular to the photometric major axis (PA ~ 45° in [O iii] and 
PA ~ 60° in Ha) and passing through the dynamic centre. Within 
our single-Gaussian model of emission profiles, the spatial re- 
gion of asymmetric lines is represented by the large-FWHM area 
of the map (Fig.[23]l. We illustrate the multi -component emission 
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profiles with the varying importance of the individual compo- 



nents in Fig. 24 



Another interesting result is that for the dust distribution. 
The map of interstellar reddening plotted in Fig.|22]shows max- 
imum extinction in the regions of emission maxima and the 
star-forming H ii ring. The obscuration of the ionisation cones 
is also asymmetric and supports the interpretation of the north- 
west cone interacting with the galactic disc, and the south-east 
cone situated above the disc. 

The stellar population ages resulting from our SSP modelling 
show a clear systematic distribution, with significant contribu- 
tions (more than a few percent) from young stars only outside 
the central ~600pc. 

NGC 5929 

NGC 5929 forms an interacting pair with NGC 5930 (e.g.. 
Bower et al.|1994| . The galaxy hosts a linear radio source with 
two lobes separated by 1" a long PA = 60° ( [Whittle et al.|1986 ' 



[Ulvestad & Wilson||1984b|l, cl ose to the galaxy major axis at 
PA = 45° ( jSchmitt et al.|1997| l. Optical emission in [Oiii] and 
Ha + [Nil] is in the form of cones co-aligned with the radio 
source, characterised by the opening angle of 80° and an ex- 
tent of 150 - 175 pc, and separated by a perpendicular dust lane 
( [Bower et al.|1994[ ISchmitt & Ki nney |1996| l. |Femiit et al.| ( [T997l l 
performed IFU observations in Ha, [Nii], and [Sii] with the 
TIGER spectrograph, and tested several models to account for 
the optical and radio emission. They place constraints on the in- 
teraction of the radio ejecta with the ambient gas, i.e., jet/cloud 
interaction, jet/ISM interaction, plasmon expansion, and the bow 
shock. None of the models were able to interpret the data satis- 
factorily. 

Our observations were completed in one spectral domain 
(MRl) only. We found differences between the major axes of 
emission in [O iii] (PA ~ 90°) and in H^ (PA ~ 75°), see Fig.|2j 
A similar discrepancy was seen in the zero-velocity curve orien- 
tation (approximately perpendicular to the emission major axes) 
in the maps of mean LOS velocities. However, the main features 
of the velocity fields, as well as their one-dimensional cuts, were 
consistent for both ions. 

Our map of stellar velocities showed strong misalignment 
with the gas velocities: the stellar kinematic axis (PA ~ 30°) was 
almost perpendicular to the symmetry axis of the gas. Moreover, 
the polarity of the stellar velocity field was reversed with respect 
to gas: in the north-east region, stars had receding motion, while 
the gas was approaching. 

The LOS velocity dispersion maps show a stripe of max- 
imum dispersion at PA ~ 160° through the optical nucleus, 
approximately perpendicular to the photometric axis of Hy6 
(Fig. 25 I. The large-FWHM region is associated with asymmet- 
ric line profiles (Fig. 25 i, with a separation of ~900kms"' be- 



tween the line components. 

The [Oiii]/Hj6 intensity ratio distribution follows the mor- 
phology of [Oiii] emission, with maxima at both the nucleus 
and off the nucleus at 2" east, 2" west, and 4" west. The stellar 
populations show a clear pattern with a growing contribution of 
young populations toward the edges of the observed field. 

5. Discussion 

Despite the advances in 3D spectroscopic techniques, optical 
datacubes obtained with the resolution that we present in this pa- 
per remain rare. We have shown the wealth of information that 



the OASIS IFU data can provide, making it possible to address 
several important problems concerning AGN, while opening a 
perspective of promising innovative, detailed modelling. The rel- 
atively large number of emission lines contained in the data not 
only allow comparisons of NLR morphologies and kinematics 
for different ionic species, but also provide diagnostics for prob- 
ing the physical conditions of the gas, based on the line intensity 
ratios. We have been able to estimate directly from the measured 
ratios both the electron density (the [S ii] doublet ratio) and the 
interstellar reddening due to dust (the Balmer line ratio). The 
determination of other characteristics such as the ionisation pa- 
rameter, metallicity, or temperature will require future modelling 
using spatially resolved diagnostic diagrams presented in this pa- 
per 

The emission data results have been compared with the stel- 
lar data, the kinematics in particular. We have also modelled 
the contributions of stellar populations of different ages and 
their variations across the observed field. This simultaneously- 
acquired, rich information allows a more complex understanding 
of the structure and conditions in each of the observed Seyfert 
galaxies. In addition to studying each of the targets individu- 
ally, we take advantage of having a uniform data set of Seyfert 2 
galaxies and search for general trends in the physical properties. 
In this section, we summarise and discuss the obtained results in 
the context of the entire sample. 

5.1. Non-axisymmetric potentials 

One of the major goals has been to map the gas velocity fields in 
detail for different emission lines, and compare them with stellar 
velocities if possible, to understand more clearly the emission 
morphologies, the dominant dynamic processes at the galactic 
centre, and eventually the origin of the NLR gas. A number of 
scenarios have been proposed for NLRs that can be tested by 
spectroscopic studies. One of them concerns departures from the 
axial symmetry of the gravitational potential, which might be 
related to the nuclear activity and the gas transport across large 
distances toward the central galactic regions (e.g., Dumas et al. 
2007] [Schinnerer et al.n 2000). On the other hand, signs of gas 
motion in non-coplanar orbits with the host galaxy may imply 
an external origin of the NLR gas, emphasising the role of merg- 
ers (e.g., [Morse et al.|| 1998} [Hunt & Malkan 2004). Sti-eaming 
motions toward the galactic centre might indicate a connection 
with AGN fuelling, and outflow motions are tested for their re 
lation to the radio jets (e.g., Bic knell et al.|199^ i, or to the pre 



Krolik 



dieted thermal wind produced by the central torus (e.g 
& Begelman 1986; Balsai-a & Ki-olik 1993) 

Only one of the gas velocity fields obtained from our ob- 
servations seems to be consistent with circular rotation in the 
central kiloparsec - NGC 449 (Mrk 1) - as far as our field-of- 
view allows us to judge. The mean LOS velocities of the ionised 
gas in seven of the eleven studied Seyfert 2 galaxies are charac- 
terised by S-shaped isovelocity contours in the central regions. 



This result supports the findings of Dumas et al. (2007 1, who re 



ported velocity contour twists to be more common in the Seyfert 
galaxies than in quiescent ones. The S-shapes are in general sig- 
natures of non-axisymmetric galactic potentials containing bars, 
warps, or nuclear spirals (Robert s et al.|1979] ), and were success- 
fully reproduced by these models in the case of several Seyfert 
galaxies e.g., by Veilleux et al.[ ( |1999a l; Schinnerer et al. ( 2000| l; 
and Emselle m et al.[ ( [2006| l. Their interpretation of the veloc- 
ity fields by motions in non-axisymmetric potentials thus pro- 
vides an alternative to the radial flows often invoked to account 
for deviations from circular velocities. Degeneracies in the inter- 
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pretation of non-circular motions were addressed analytically by 
Wong et al. ( 2004 1, who derived prescriptions of how to discrim- 
inate between bars, warps, and radial flows by means of Fourier 
analysis. 

Beside the S-shaped velocity fields, another type of a non- 
axisymmetric velocity pattern has been represented in our data 
in the case of NGC 262 (Mrk 348): two striking spots in the mean 
velocity map, suggesting a rotating disc/ring (~ 300 pc in radius) 
that is non-coplanar with the host-galaxy disc. Concerning out- 
flows, we cannot provide direct evidence, even though at least 
two of the objects in our sample (NGC 4388, NGC 2992) have 
had an outflow component confirmed by several authors (e.g., 
[Marquez et aT][TW8l |AUen et al.|[lW 



Veilleux et al. 1999b 



Veilleux et al.|2001[ ). We are partially limited by the small FOV, 



but a clearer understanding is expected from our future detailed 
kinematic modelling. 

Except for one case (NGC 4388), the non-circular motions of 
gas do not have evident counterparts in the stellar velocity fields. 
This is unsurprising since gas, even in the absence of an AGN- 
driven outflow, is in general not expected to follow the same 
orbits as stars, due to its dissipative nature resulting in lower ve- 
locity dispersion, its higher responsiveness to non-axisymmetric 
perturbations, possible ring formation, and a phase shift with re- 
spect to axes of stellar orbital families. A counter-rotation or tilt 
with respect to the plane of the stellar disc are also plausible in 
the case of gas having an external origin (e.g., as a result of a 
minor merger). Observed misalignments between the kinematic 
axes of the stellar and gaseous components in active and non 
active galaxies were discussed in detail by [Dumas et al.| ( [2007 



who concluded that departures are both more frequent and sig- 
nificant in galaxies hosting an AGN. They also found that the 
closely aligned axes and a lack of isovelocity contour twists in 
gas were correlated with low-accretion rate AGN. Their results 
suggested a close relation between the AGN and the dynamics 
of the host. 



5.2. Nuclear rings 

Further type of evidence of a non-axisymmetric potential is the 
presence of nuclear rings. A small fraction of observed rings 



may be due to collisions or mergers of galaxies (e.g., Knapen 
et al. 2004[ l, or to accretion of intergalactic gas, although the vast 



majority are probably resonance phenomena ( Buta & Combes| 
19^ . Nuclear rings have been identified in approximately 20% 
of spiral galaxies ( | Knapen 2005 i. Our sample of Seyfert 2 galax- 
ies possesses at least one object with a nuclear ring, reported 
previously ( [Wilson et al. 1993| l and confirmed in our data: 
NGC 5728. We detect the ring in the surface brightness maps 
of the Balmer lines and low-ionisation lines of metals, while 
it is missing in [Oiii]. We have found it to be associated with 
high interstellar extinction, which suggests that it contains large 
amounts of dust. 

Evidence of other rings might originate in the kinematic 
measurements, as was the case for NGC 5252 in [Morse et al.| 
( [199 8), and will be tested in future modelling. The strongest 
candidate to contain a ring is NGC 262 (Mrk 348) with the strik- 
ing velocity pattern describ ed in Sect.[5.1[ n ot reported so far 
elsewhere. As suggested by Mulchaey et al. ( 1996a| l, NGC 262 
possibly hosts another ring at ~ 1 kpc. Marginal evidence of its 
detection was also found in our data, in the mean LOS velocity 
map of [Oiii], which displays systematic velocity pattern at the 
edges of our FOV that is qualitatively consistent with rotation in 
an opposite sense to that of the nuclear ring. 



5.3. NLR morphologies, cones 

The morphologies that we have obtained from analysis of dif- 
ferent individual lines are not always identical or even aligned 
along the same axis. The misalignments between emission in 
[Oiii] and low-ionisation fines (e.g., NGC 449, NGC 5929) 
provide evidence of additional dynamic processes affecting the 
high-ionisation components of the gas, such as the radio jets that 
accelerate and/or ionise the gas. 

Two of the galaxies (Mrk 1066, NGC 5728) show prominent 
regions detectable only in hydrogen and low-ionisation lines, 
which are virtually missing in [O in] : a nuclear ring in NGC 5728 
(reported before e.g., by Wil son et al.j 1993"] l, and an off-centre 
region of diameter ~500pc in Mrk 1066. The regions are char- 
acterised by a ratio [O iii]/Hj8 markedly lower than in the remain- 
der of the observed field, and have therefore probably not been 
ionised by the AGN radiation. They are also associated with high 
interstellar reddening, which indicates a high amount of dust and 
possibly star formation, and hence an ionisation by young mas- 
sive stars. Marginal support of this interpretation comes from 
the maps of stellar ages, which show relatively large contribu- 
tions from young stars in the specified locations compared to 
the neighbouring regions. Ionisation by sources other than the 
active nucleus thus definitely play a role in the observed NLR 
morphology. 

Direcdy related to the Unified Model of AGN (and the ex- 
istence of the nuclear torus), the major puzzles about NLRs 
are their extent, the existence of ionisation cones, and the in- 
trinsic geometry of the gas distribution. Significant evidence of 



cones was found in ground-based and HST images (e.g., [Wilson 
et al. 1993; Falcke et al. 1998). The NLR structure parameters 
were analysed statistically for large samples by Mulchaey et al. 
( [1996a| l; and |Schmitt et al.| ( |2003 a), who both concluded fliat a 
surprisingly low number (with respect to the predictions of the 
Unified Model) of NLRs in Seyfert 2 nuclei have a major axis 
that is significantly longer than their minor axis. We reached the 
same conclusion, with the caveat that our results are influenced 
by both the small FOV and spatial resolution, mostly determined 
by the seeing conditions (~1")- 

[Mulchaey et al. ( 1996b) made predictions of NLR extents 
and morphologies for different observer viewing angles, and also 
for different intrinsic gas distributions ionised by the conical ra- 
diation field from the AGN. They considered two basic geome- 
tries of the gaseous region, a sphere and a disc, and tested their 
model predictions with data from a large imaging survey. In this 
context, our kinematic measurements clearly show that the ob- 
served NLRs are characterised by mostly rotational (although 
probably non-circular) motions, and are thus non-spherical. 

5.4. Alignment with radio emission 

Where detected, the elongated optical emission structures have 
been found to be aligned with radio emission data published in 
the literature, to within the measurement errors. This alignment 
infers an alignment between the torus and accretion-disc axes 
and/or an influence of the radio jet on the NLR gas dynam- 
ics and ionisation ( [Schmitt et al.[[2003b] l. Our result confirms 
previous studies of the alignment between the two wavelength 
bands, demonstra ted e.g., by|Wilson & Tsvetanov ( l994j; Nag^] 



Schmitt et al.[ ( [2063b| ). 



& Wilson[([T999l);[Falcke et al.[ ( [1998| l; |Capetti et a 



( 1996 1; and 



As discussed above, NGC 449 and NGC 5929 of our Seyfert 
sample have different orientations of photometric axes in the 
high-ionisation and low-ionisation emission lines. The radio-jet 
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axis is consistent with [O m] emission in NGC 449, and with Hy6 
emission in NGC 5929. Conclusions are ambiguous for NGC 
2273, which exhibits different shapes in its radio continuum 
emission at 20 cm and 6 cm wavelengths ( jUlvestad & Wilson] 
|1984b| ). The Hf) emission appears to agree with the former and 
[O III] with the latter. 



5.5. Emission-line profiles 

We detect prominent splitting of emission lines in eight of the 
observed NLRs, usually limited to only part of the FOV, iden- 
tified as a large-FWHM region in the maps produced by our 
single-Gaussian line-profile modelling. In seven of the galax- 
ies the region is of elongated shape approximately perpendicular 
to the NLR major axis and passes through the dynamic centre. 
Multi-component emission-line profiles have been reported in 
numerous NLR spectroscopic studies, but the integral-field spec- 
troscopy permits us to map in detail the spatial variations in each 
component and to use the information in modelling the NLR 
structure and kinematics. The origin of the line splitting found in 
our data is unclear and will be addressed in future modelling of 
the observed objects. The most common interpretations include 
superposition of several kinematic components (e.g., Schommer 
|et al. 1988 ), or expanding outflows (e.g., ICapetti et al. I999[ l, 
rapidly rotating circumnuclear discs (e.g., Rubin et al. 19 97[), and 
also bar-induced non-circular motions (e.g.,|Veilleux et al. 2001 



Bureau & Athanassoula|1999||Kuijken & Merrifield|1995| . 



5.6. Electron densities 

We have estimated electron densities from the [S ii] doublet in 
six galaxies of our sample for which data are available. In all 
cases, a non-uniform density distribution with a central peak has 
been found, the peak values ranging between ~ 600 cm"^ and 
~1500cm"^, and declining to ~ 200 cm"^ outside the central ~ 
500 pc. The derived inhomogeneous gas distributions may have 
a strong effect on the estimates of global quantities, such as the 
NLR mass and the Stromgren radius. 

The distributions are asymmetric and the aforementioned 
values are thus direction-dependent, although the overall results 
are consistent with long-slit observations (e.g., |Bennert et al.| 
2006a Fraquelli et al. 2003] ). More precise modelling is nec- 
essary to decide which of the values in the outer parts of the 
FOV (which are often relatively high) estimate reliably the in- 
trinsic densities and which are an effect of measurement er- 
rors, to which the derived density values are highly sensitive. 
We also note that both the density values and the radial pro- 
files are influenced by the method of derivation. The [S ii] ratio 
is density-sensitive only within a limited range of values, and 
provides a density measure only in regions where these sulphur 
ions are present. The effects of observational errors are also non- 
negligible. We have also derived the densities based on the as- 
sumption of an electron temperature of T - lOOOOK, and dif- 
ferent temperatures in different regions may therefore modify 
the absolute values and the spatial distribution of the computed 
densities (see e.g., ]Osterbrock|19891 l. 



5. 7. Interstellar reddening 

The dust distribution may contribute significantly to the recon- 
struction of the galactic structure and the orientation of the ob- 
served nuclear emission regions with respect to the host galaxy, 
such as emission cones above or underneath the galactic disc. 



We have been able to estimate the dust distribution from the ra- 
tio of Balmer lines in 7 objects. The maximum extinction is usu- 
ally connected with regions of maximum emission. Distinct off- 
nuclear regions of high extinction have been found in NGC 262 
(Mrk 348) and NGC 5728, consistent with low-ionisation re- 
gions such as the star-forming ring in NGC 5728. The dust dis- 
tribution in NGC 5728 is an excellent example of asymmetry: 
the obscuration is much higher in the north-west emission cone 
than in the south-eas t part (see Fig.]22| , which is consistent with 
the interpretation of ]Schommer et aL| ( 1988) 1, and Wilson et al. 
( 1993) 1, who predict the north-east cone to be interacting with the 
galactic disc. 

The values of E(B - V) that we have derived were based on 
the assumption that the intrinsic Ha/HyS ratio value equals 3. 1, as 
is standard for AGN-excited gas due to hard radiation causing a 
relatively large transition zone ( Osterbrockj 1 989) 1. However, not 
all of the observed spectra necessarily originate in AGN-excited 
locations with an overabundance of high-energy photons com- 
pared with stellar spectra, and the intrinsic value of 2.86 would 
have been more appropriate. Therefore, reddening is underesti- 
mated in these locations. 

The E(B - V) that we have obtained is strongly dependent on 
the stellar population model applied to correct the Balmer emis- 
sion lines for underlying absorption. The computed E(B-V) val- 
ues are significantly higher than those published e.g., by Bennert 
etal., (2006a); Raimann et al. ( 2003) ; Fraquem"etaLl ( ]20031 l; and 
Ho et al.l ( 11997) l. The values that are common in the literature are 



of below one magnitude, in many cases determined by methods 
that differ from that of the Balmer decrement, and are thus in- 
dependent of the stellar absorption models. We have tested the 
influence of the SSP models on the spatial distribution of dust, 
and on the E{B - V) values. We have found that the spatial distri- 
butions remain unchanged qualitatively after applying the "min- 
imum absorption correction" and using an old SSP as a stellar 
template. In contrast, the distribution differed when we did not 
correct for absorption and subtracted only the smooth contin- 
uum. We considered the no-absoiption case to be non-realistic, 
and therefore conclude that our maps of E{B - V) reflect the true 
character of the dust distribution. At the same time, the prob- 
lem of the high reddening values was not solved by applying a 
different stellar model. 

5.8. Diagnostic line ratios 

The maps of diagnostic line ratios trace the spatial distribution 
of the gas physical characteristics and the properties of the ioni- 
sation source. The [O iii]/Hy6 ratio provides the most straightfor- 
ward interpretation by tracing the gaseous media ionised by dif- 
ferent sources projected along the line of sight. An outstanding 
example is NGC 4388 in ]Veilleux et"aL] ( ]1999b| l, where [Oiii]/Hje 
helps to separate between the gas in the galactic disc and the ex- 
traplanar material. An impact of the ionisation maps has been 
demonstrated in the case of NGC 5728, where the higher ioni- 
sation in one of the lobes has been used to advocate the inflow 
interpretation of the observed velocity field due to the position 
of the ionisation cones in relation to the galactic disc ( [Schomrrier 
ietal..l988 ). 

In the objects that we have studied, maximum ionisation oc- 
curs in the vicinity of the galactic nucleus, and is consistent 
with a model of ionisation by AGN radiation. The isocontours 
traced by the [Oiii]/Hj6 are not identical to those of the [Oiii] 
emission, and in approximately half of the sample the symme- 
try changes, which implies a contamination by other emission 
sources. The ionisation maps in general show more clearly de- 
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fined axes than the [O m] emission, and have a more conical form 
than the [Oin] maps, as far as our spatial resolution allows us 
to perceive. As already discussed in previous subsections, the 
ionisation maps have helped us to identify a ring in NGC 5728 
and probably a star-forming region in Mrk 1066. We have also 
confirmed high ionisation in the southern lobe of NGC 5728, re- 
ported by |Schommer et al.| ( [T988 ). The [Oil/Ha, [Nii]/Hq', and 
[S ii]/Ha ratios have distributions that are ofF-centre in most of 
the observed cases. 

We have constructed a set of three spatially resolved diag- 
nostic diagrams for 7 of the 1 1 observed Seyfert 2s (for Mrk 34 
only two have been possible). The distributions of points in 
the diagrams show considerable differences between the NLRs, 
reflecting important differences in their inner structures. We 



have included several theoretical curves (Kewley et al. 



Kauffmann et al. 2003 ; Stasihska et al. 2006 Kewley et al. 



2001 



2om 



classifying emission-line objects based on their photoionisation 
source. However, the division inside the plots is only informa- 
tive, since the models correspond to spectra integrated spatially 
over the nebula and are not adapted to the spatially resolved dia- 
grams. A numerical ionisation modelling is necessary to achieve 
a correct interpretation. A strong Seyfert regime across the entire 
observed range of projected radial distances is detected in Mrk 
34 and NGC 262 (Mrk 348). The opposite is true for Mrk 622, 
where all of the points lie below the Seyfert/LINER limit (which 
may be a consequence of the applied SSP models). The distribu- 
tions of points are much less compact for Mrk 1066, NGC 449, 
NGC 2992, and NGC 5728, which have regions in the Seyfert 
regime as well as regions consistent with non-Seyfert ionisa- 
tion. Numerical modelling similar to that carried out by Bennert 



et al.| p006a| ) will be necessary to determine the true edges of 



the NLRs. 



5.9. Stellar populations 

Stellar population modelling performed by using synthetic 
young and old SSPs has revealed varying contributions from 
young stars and their spatial distribution. The occurrence of 
young stars confirms results achieved in previous studies by sev- 
eral independent methods: wind fines from O and B stars in 
ultraviolet (UV) spectra (e.g., IHeckman et al. [1997 Gonzalez 
[Delgado et al. |1998| l, and absorption fines of fielium and of 
fiydrogen in fiigfi-order B aimer series from O, B, and Wolf- 
Rayet stars (e.g., |Gonzglez Delgado et al.ll998 l|2001|). Statistical 
studies of Seyfert galaxies (e.g., Storc hi-Berg mann et al-POOO" 



Joguet et al.||2001 Raimann et al.||200 3 ) reported signatures of 
young populations in 30 - 50% of the observed Seyfert 2 galax- 
ies, and significantly higher contributions (~ 20% in terms of 
flux in the nuclei) than in non-active galaxies. The detection of 
recent starburst in the nuclear regions of active galaxies may be 
closely related to the interplay between the host galaxy and the 
AGN, and to the triggering of the AGN . The highest contribu- 
tions from young stars are present in AGN with high IR luminos- 
ity, which are interacting and have distorted inner morphologies 
( Cid Femandes et al.|2001 Storchi-Bergmann et al.|200T i. 

We have compared the relative mass fractions of young stars 
obtained in our study with the results of |Cid Fernandes et al. 
( |2004| l. The comparison has been completed for spatially unre- 
solved data, where we have simulated a 1" aperture to comply 
with |Cid Fernandes et al.| ( |2004| l (see Table[4]l. The mass frac- 
tions that we obtain for the young population are of the same or- 
der as the combined fractions of young (< 10^ My r) and interme- 
diate (10^ Myr - 1 .4 Gyr) populations inferred by Cid Femandes 
^taL] ( [2004l i. 



Table 4. Weights of stellar populations obtained by modelling of 
OASIS data, and compared tol Cid Femandes et al.| ( |2004) , on aperture 
~l" . The symbols fiY,fii,fio denote mass fractions of young, interme- 
diate and old populations, respectively. Analogously for /Jyounj, Hou in 
our OASIS data. 



Galaxy 


OASIS data 


Cid Femandes et al. 1 2004 1 


A'young A'old 




NGC 2992 
NGC 3081 
NGC 4388 
NGC 5728 


0.09 0.91 
0.008 0.992 
0.02 0.98 
0.06 0.94 


0.0000 0.0260 0.9740 
0.0000 0.0000 1.0 
0.0002 0.0000 0.9998 
0.0002 0.0699 0.9299 



At least half of our Seyfert 2 sample exhibit clear spatial gra- 
dients in stellar population ages. In four galaxies (NGC 3081, 
NGC 4388, NGC 5728, NGC 5929), old stars have been found 
to be dominant close to the dynamic centre or the location of 
peak emission, and a gradual increase in the contribution from 
young stars was measured outwards, which contradictsiRaimann 



et al. (2003) who found an opposite trend. On the other hand, 
NGC 2992 contains its youngest populations close to the dust 
lane, and NGC 2273 has a maximal contribution from the young 
SSP north of the dynamic centre, which gradually decreases out- 
wards. The remaining five galaxies (NGC 449, Mrk 34, NGC 
262, Mrk 622, Mrk 1066) do not have high-quality stellar data, 
but the youngest populations have been fitted visibly in the 
central-most regions. 

6. Summary and conclusions 

We have analysed a set of 11 nearby Seyfert galaxies of type 
2, observed with the integral-field spectrograph OASIS, which 
mapped the central ~ 1 kpc with medium spatial and spectral 
resolution. The central goal of our study was to obtain - 
from optical emission-line spectra - spatially resolved spectral 
diagnostics, and kinematical maps of the ionised gas. We have 
constructed, where possible: 



2D maps of: surface brightness in individual lines, line 
ratios, velocity fields (mean LOS velocity and velocity 
dispersion), electron density, interstellar reddening, stellar 
mean LOS velocities, and relative mass fractions of young 
stellar populations; 

Spectral-diagnostic diagrams for line ratios sensitive to the 
source of ionisation (Seyfert, LINER and starburst regimes). 



The 3D technique that we have applied provides a wealth 
of information that will be investigated in depth in forthcoming 
papers. The combined advantage of the relatively high spatial 
resolution and the length of the spectral range produce unprece- 
dented details about the morphology, and both the kinematic and 
ionisation structure of the ionised gas of NLRs. We draw the fol- 
lowing conclusions from the present study: 

- Eighty percent of the observed galaxies exhibit S-shaped 
isovelocity contours in the mean LOS velocity maps of the 
ionised NLR gas. This suggests that non-axisymmetric po- 
tentials are probably common within the studied sample. 

- We have detected a possible nuclear ring or a radial flow of 
gas in NGC 262 (Mrk 348), which was not reported before. 



Stoklasova et al.: OASIS integral-field spectroscopy of the central kpc in 11 Seyfert2 galaxies 



15 



It has been identified based on the mean LOS velocity maps 
derived from the observed emission lines. 

- Electron densities usually reach their maximum values of ~ 
1000 cm"-' at the galactic nucleus, and decrease outwards. 

- Dust is distributed unevenly in the central regions and will 
have to be included in a future detailed study of conditions 
in NLRs. 

- Emission-line splitting has been identified in 8 out of 1 1 ob- 
jects, confined to a clearly defined spatial region at the galac- 
tic centre, and in most cases elongated perpendicular to the 
major axis of emission. 

- Our spatially resolved diagnostic diagrams show radial evo- 
lution from the Seyfert-like regime in the central regions to 
other regimes at larger distances. 

- We find spatial gradients in stellar population ages, and con- 
firm the presence of young stars. 
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velocity dispersion (FWHM, in kms line intensity ratios 
[Oiii]/H/3, [Nii]/Ha, [Sii]/Ha, [Oi]/Ua, reddening E(B - V), 
and electron density; 

b) Spatially resolved spectral-diagnostic diagrams; 

c) Selected emission-line ratios as functions of projected 
galactocentric distance; 

d) ID cuts through mean LOS velocity and electron density 
along specific position angles; 

e) Spectra - for selected individual spaxels, and/or integrated 
over circular apertures or the entire FOV - and their fits. 

In the maps, north is up, east to the left. The black contours 
represent the surface brightness in the corresponding emission 
line for each of the derived quantities: in the case of mapped line 
ratios, the contours correspond to the forbidden line included in 
the ratio; for the electron density, they represent the sum of the 
two sulphur lines; and for the maps of interstellar reddening, they 
depict the surface brightness in Ha. Any discrepancies between 
contours of the same species in different maps result from their 
re-computation for the new grid, where points of low S/N in the 
other relevant line(s) were removed. 

In the one-dimensional cuts, the projected radial-distance co- 
ordinates have been defined as increasing from east to west. 
Error bars in diagnostic diagrams and the emission-line ratio 
scatter plots have not been assigned to individual points, for 
clarity. Instead, median and maximum errors have been plot- 
ted for each of the defined spatial regions, in the correspond- 
ing colours. The spatially resolved spectral-diagnostic diagrams 
compare several classification schemes commonly applied to 
emission objects (originally with spatially non-resolved spec- 
tra). The two parallel lines for [O i]flia correspond to the values 
provided by .Ho et al.| ( [T997l ), log([Oi]/Ha') > -1.1 for Seyfert 
galaxies, and log([0 iJ/Hor) > -0.77 for LINERs. 



Caption to Figs. [2p5| 

Figures [2p5] contain (where available and relevant to the 
description of results and discussion, depending on individual 
galaxies): 

a) 2D maps of surface brightness in emission lines of 
Ha, H/3, [Om], [Ni], [Oi], [Nii], and [Sii] (in units of 
arcsec ^m ^), mean LOS velocity (in kms '), 
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Fig. 2. Mrk 34. See detailed caption before this graphic section. Surface brightness is in 10 J s ' arcsec ^ m ^, extinction in mag 
in all figures. 
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Fig. 3. Mrk 34. Mean LOS velocities and FWHM are in unites of km s ' in all subsequent figures. 
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Fig. 4. Mrk 622. 



Mrk 622 - Ratio [5 II] 6717,6731/ Ha 




Diagnostic diagram |0 lll]/Hbvs. |0 l]/Ha 



++H- I- < T <2- 
r < l- 

KMliY^dl.iaiL 

hki ^ri. 

hki ^ri. ^cvfirt 

hki ^ri. 







1-3-1 


i-i-i 




-i.O -1.3 -1.0 -03 

log f [0 I] eSDD/ H_alpha J 



20 Stoklasova et al.: OASIS integral-field spectroscopy of the central kpc in 1 1 Seyfert 2 galaxies 

Mrk 622 - Velocity in [0 Nil r^lrk 622 - Velocity in Hb Mrk 622 - Velocity in Ha 






I 


—!i i- 


* i i 


1 








arcsec 




440 












Mrk622 


- FWHM of [0 III 






-!s * 1- 



-160 



J . 

u 
1/1 

u 0- 




Mrk 622 - FWHM of Hb 





S3 -160 



Mrk 622 - FWHM of Ha 



S3 




900 



900 



600 



Mrk 622 - Electron Density 




100 



10000 



Mrk 522 - Vslocitv PA-Q.Q 



SSP mass Young/Total, S/N=50 




Mrk 62 2 - ElEctran d En sit/ PA-Q.Q 




Mrk 522 - Diagnasbc ratia IN II] 6533 / Ha 




Mrk 522 - Diagnostii ratio [S II] 5717,5731 / Ha 



+++I- r < L- 

-- hfci Jt J. L3J7 



ii 



R larcsEc] 



f'^fe — ti ta — ti — ta — lis i4 -tr 

R larcsEc] 



Stellar velocities S/N=50 




Mrk 622 - Diagnastic raba |0 III] 5QQ7 / Hb 




if 



'■]] 



1 i i 



R larcsEc] 

Mrk 62 2 - Diagnastic ratia |0 I] 63 QQ / Ha 



i- 



.(.-i ; — 

* J- <r ^ 1- 



1- i-t* 

.^.t4... 







u 


11 



isnfaar 

H II 



llil iLi [.il I.S lil lit 

R I arc sec] 



Fig. 5. Mrk 622. 
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Fig. 6. Mrk 622. 
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Fig. 7. Mrk 1066. 
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Fig. 8. Mrk 1066. 
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Mrk 1066 - integrated MRl spectrum 
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Fig. 9. Mrk 1066. 
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Fig. 10. NGC 262 (Mrk 348). 
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Fig. 11. NGC 262 (Mrk 348). 
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Fig. 13. NGC 449 (Ivlrk 1). 
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Fig. 14. NGC449 (Mrk 1). 
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Mrk 1 - integrated MRl spectrum 
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Fig. 15. NGC449 (Mrk 1). 
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Fig. 16. NGC 2273. 
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Fig. 17. NGC 2992. 
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Fig. 18. NGC 2992. 
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NGC 2992 - integrated MRl spectrum 
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NGC 2992 - MRl, asymmetric lines 
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NGC 30S1 - Intensity in [0 III] 5007 
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Fig. 20. NGC 3081. 
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Fig. 21. NGC 4388. 
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Fig. 22. NGC 5728. 
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Fig.23. NGC 5728. 
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NGC 5728 - integrated MRl spectrum NGC 5728 - integrated MRl spectrum 
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Fig. 24. NGC 5728. 
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Fig. 25. NGC 5929. 



